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The High-Density Symmetry Energy and Direct Urea 

Andrew W. Steiner^ 
^ Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545, USA 

The symmetry energy of nucleonic matter is usually assumed to be quadratic in the isospin density. 
While this may be justified at sub-saturation densities, there is no need to enforce this restriction at 
super-saturation densities. The presence of a quartic term can strongly modify the critical density 
for the direct Urea process which leads to faster cooling of neutron stars. Neutron star cooling 
predictions which lie below the observational data can, for some equations of state, be repaired with 
a quartic term which effectively turns off the direct Urea process. 

PACS numbers: 97.60.Jd, 26.60.+C, 21.65.+f 



O 
O 

(N 



oo 



> 

O 
O 

o 
\o 
o 



X 



I. INTRODUCTION 

The equation of state of nucleonic matter is often de- 
composed into an isospin-symmetric and an isospin asym- 
metric part 



E{n, 6) = Enuc{n, 6 = 0) + S^S{n) + 6^Q{n) + 



(1) 



where n = n„-|-np,J=l — 2np/ n, i?nuc is the energy per 
baryon of nuclear matter and S and Q are arbitrary func- 
tions of density. Most often, Q is sufficiently small so that 
it can be ignored, and then -Bsym — S{n) is the density- 
dependent and isospin-independent symmetry energy. 

Without the Coulomb interaction, the near isospin 
symmetry of QCD predicts that pure neutron matter and 
pure proton matter should have the same energy density. 
Terms in Eq. ^ with odd-powers in delta are thus for- 
bidden, assuming the 6 expansion is analytic. Near the 
nuclear saturation density, it has been suggested that the 
truncation to order 6'^ for the potential energy part of the 
equation of state is a very good approximation for den- 
sities below or near the nuclear saturation density, 0.16 
fm-3 H S S S P, "6^. In this article, the effect of the 
quartic term in Eq.^at higher density is examined, where 
no constraints on the quartic terms are known because 
the uncertainty contained in many-body effects. 

The major effect of a "quartic term" (i.e. a term that 
is of fourth-order in the isospin deviation, 6) is to change 
the ratio of protons to neutrons at higher densities. All 
other things remaining equal, the effect on neutron star 
masses and radii from a quartic term is typically small. 
Howeverj_the critical density for the direct Urea cooling 
process 



' p + e + Ve and p + e 



(2) 



is very sensitive to the number of protons present in dense 
matter. This process can only proceed when the condi- 
tions for momentum and energy conservation can be si- 
multaneously satisfied, when the proton fraction is larger 
than about 10 percent. The presence of a quartic term 
in Eq. ^ can either move the critical density for the di- 
rect Urea process to much lower densities, or to densities 
larger than the central density of the maximum mass 
neutron star. 



Neutron stars cool much more quickly when the di- 
rect Urea process is allowed 0, . Thus understanding 
the possible magnitude of quartic terms is essential in 
describing the cooling of neutron stars. 



II. ORIGINS OF QUARTIC TERMS 

For the present discussion, it is useful to separate the 
equation of state into contributions to the potential and 
kinetic energy. While this demarcation is not necessarily 
unique, we use the term "kinetic energy part" to refer to 
the Fermi gas portion (including the interactions present 
in the nucleon effective mass) of the equation of state 
(whether relativistic or non-relativistic) and "potential 
energy part" to refer to the remainder. The kinetic 
energy part of the symmetry energy naturally contains 
quartic terms in all equations of state. The potential 
energy part of the equation of state is the primary con- 
sideration here, but it should be noted that interactions 
can affect the size of quartic terms in the kinetic energy 
through their effect on the nucleon effective masses. 

In relativistic mean-field models [TollllllTllTllTi| the 
many-body interactions have a non-trivial isospin depen- 
dence even in their most trivial form. The linear Walecka 
model, for example, already has quartic terms (albeit 
small) in the potential part of the energy density. These 
quartic terms can be much stronger when three- and 
higher-body contributions are included. An example is a 
term proportional to ap^ where a is the scalar-isoscalar 
meson and p is the vector-isovector meson. In the context 
of a relativistic point-coupling model 0, IT^ , such a 
term would correspond to a three-body nucleon-nucleon 
interaction of the form 



(3) 



In Skyrme-like models, no quartic term is present in 
the potential energy part of the equation of state. The 
two-body interaction present in the Skyrme model results 
in two powers of density, and the three-body interaction 
is treated as a density-dependent two-body interaction 
where the additional density dependence is assumed to 
be exactly isospin symmetric. 
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Modern microscopic-macroscopic and Hartree- 

Fock-based JJj models of nuclear masses contain explicit 
non-quadratic terms in the form of the phenomenologi- 
cal term called the Wigner energy. This term is often 
thought to originate in isovector pairing pft. .2L ..22. .231 
0,|2^. Pairing is thought to be primarily a surface effect 
and vanishes in the infinite baryon number limit. Thus it 
does not contribute to the equation of state of infinite nu- 
cleonic matter here. However, if the Wigner energy was 
shown to have a different origin, or if the pairing energy 
was shown to be partially a volume effect (rather than 
a surface effect), then it would contribute to the isospin 
dependence of the symmetry energy. Neutron star mat- 
ter, however, is typically sufficiently neutron-rich that 
isovector pairing is not a strong contribution. 



III. CONVENIENT PARAMETERIZATION OF 
THE QUARTIC DEPENDENCE 

A standard alternative to Eq. ^ is to define the sym- 
metry energy using the second derivative with respect to 
the isospin asymmetry 



£^sym(n,(5) = - 



f (fE{n,6) 



= S{n) + 6Q{n)d^ , (4) 



where e is the energy density. In this work, Eq. ^ will 
be taken to be the definition of the symmetry energy. If 
there are no quartic terms in Eq. ^ then iJsym — S{n) is 
independent of the isospin asymmetry, S. In the presence 
of quartic terms, this definition of the symmetry energy 
in Eq. ^ is more closely connected to the experimental 
observables based on nuclei which like close to the N = Z 
line than the alternative of defining the symmetry energy 
as the energy difference between pure neutron matter and 
isospin-symmetric nuclear matter. 

A convenient way to parameterize the addition of a 
quartic term is to define rj as 



ri{n) 



E{n,l) - E{n,l/2) 4S{n) + 5Q{n) 



3[E{n,l/2) - E{n,0)] AS{n) + Q{n) 



(5) 



This parameter is unity when the quartic term is zero 
and greater or less than one depending on whether the 
symmetry energy is increased or decreased by the quartic 
term. This parameter is more convenient for numerical 
work than using the generalization of Eq. 0] , 



E. 



sym , 4 



(n) 



1 rf4£;(n,(5) 

24 dS^ 



Q{n), 



(6) 



since numerical fourth derivatives can be difficult to com- 
pute accurately. The value of rj is restricted by ensuring 
that the derivative of the energy per baryon as a function 
of delta does not vanish at any point other than aX 6 = 
corresponding to isospin-symmetric nuclear matter. This 
guarantees that, at fixed density, nuclear matter is always 
the most energetically favored configuration and neutron 



matter is always the most costly configuration (when ig- 
noring the Coulomb interaction and possible constraints 
from /3-equilibrium) . The proper hmits on -q, 3/7 < rj < 5 
are easily obtained from the definition above. Sometimes 
it is useful to distinguish the value of rj as obtained from 
the "kinetic" part of the EOS, T^kin, that obtained from 
the "potential" part of the EOS, jypot, and that obtained 
from the full EOS, jytot- Note that ?7tot is not trivially 
related to the individual contributions Tykin and ?7pot- 

Enforcing a particular density-dependence of rj on the 
potential part of a given equation of state is straightfor- 
ward. Defining T as 



[~3xp + 19x1 - ^'^4 + ^^4+ 



1 + 377 

7xpf] - 23Xpr] + 32x^7] — IGx^i]] , (7) 



where Xp = Up/n, one can construct a new equation of 
state from any equation of state of neutron and nuclear 
matter using 

Epot{n, S) = T£'pot,nuc("-) + (1 - ^)i^pot,nout(") • (8) 

(Note that the limit rj —> 1 gives the correct expression 
for a purely quadratic symmetry energy, 1 — (1 — 2xp)^ .) 
Only the potential energy part of the equation of state 
has been modified, and the kinetic energy part remains 
unchanged. The quantity T is density dependent if rj is 
density dependent, as is typical for many equations of 
state. The corresponding potential parts of the neutron 
and proton chemical potentials are given by 



Mpot,n — ^/^pot.nuc.n ,ncut,n ~r 

,nuc,n J 



, {Xp 1) , 

x'^P ) VP'pot,ncut,n 



+ (1 - T)^pot ,nout,p ~r 
T'{Xp)-^ (Mpot ,ncut,p /^pot,nuc,p ) (9) 

The effective masses are not modified here since they are 
in the kinetic energy part of the EOS. 

The behavior of the energy per baryon as a function of 
5 at saturation density in the equation of state (EOS) of 
Akmal, et. al (APR). (26j] constructed as above to have 
varying values of 77pot is given in Fig. ^ The deviation is 
not large, at most a couple MeV near saturation density 
for the most extreme values of rj. 



IV. THE VALUE OF r/ FOR MODERN 
EQUATIONS OF STATE 

In addition to APR from above, several relativistic 
mean-field (RMF) models are used in this work, all of 
which are constrained to have the following properties at 
saturation 
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FIG. 1: The energy per baryon of the APR EOS at saturation 
density as a function of the proton fraction, Xp, for the trivial 
and limiting values of r^pot- 

K = 230 MeV 
Esym = 34 MeV 
M7M = 0.8 (10) 

where np is the saturation density, E/A is the binding en- 
ergy, K is the compressibility, Esym is the symmetry en- 
ergy at the saturation density, and M*/M is the reduced 
Dirac effective mass. In addition to the standard RMF 
model with no non-linear couplings except for the cubic 
and quartic self- interactions among the scalar mesons, 
two RMF models are constructed with small ( "RMFlo" ) 
and large ( "RMFhi" ) values of ?7pot in order to demon- 
strate the variations that are possible within this formal- 
ism. These two models do not necessarily represent the 
most extreme cases, as a full optimization was not per- 
formed. The couplings are given in Table 1 and utilize the 
Lagrangian and notation described in |23|. AH of these 
couplings are within a factor of two of the constraints 
required by naturalness 13, 28] . 

The Skyrme ^ model "SLy230a" ^30l] is used for com- 
parison. This model matches the binding energies and 
charge radii of several ground-state nuclei, the known 
properties of saturated nuclear matter, and results in rea- 
sonable neutron star properties. 

The momentum dependence of the nucleon optical po- 
tential is essential for transport simulations of heavy-ion 
collisions. The MDI (momentum-dependent interaction) 
EOS 31], which has successfully described aspects of 
heavy-ion collisions at intermediate energies, is included 
here. In addition to providing for reasonable properties 
of saturated nuclear matter, the isoscalar potential of the 
MDI EOS coincides with predictions from the variational 
many-body theory using inputs constrained by nucleon- 
nucleon scattering data 0, l32| , and the isovector poten- 
tial agrees with the momentum dependence of the Lane 



TABLE I: The couplings for the RMF, RMFlo, and RMFhi 
models. The units are all given so that the Lagrangian has 
units MeV* when the meson fields are in MeV and the nucleon 
fields in MeV''''^. Omitted couplings are equal to zero. 



RMF 


5a 


7.721 


9u, 


7.955 


9p 


8.608 


K 


21.40 


A 


-6.227 xlO"* 






RMFlo 




7.692 


9u> 


7.998 


9p 


14.00 


K 


18.9 


A 


5.188 xlO"^ 


c 


0.03801 




1.499 




69.6 




1.052 


bi 


0.05251 


as 


8.446 xlO"^ 




-2.063 xlO"^ 




-8.085 xlO"*^ 






RMFhi 


5t 


7.677 


5" 


8.021 


9p 


7.551 


K 


17.50 


A 


8.361 xlO-3 


c 


0.05880 




0.02870 


Ol 


-35.30 




-0.66605 


bi 


0.2719 


as 


1.281 xlO"^ 


a4. 


2.062 xlO~^ 


&2 


1.127 xlO"'^ 







potential extracted from (p,n) charge exchange reactions 
up to about 45 MeV [s^ 13^ . We use the formulation as 
described in Ref. [sij with the symmetry energy param- 
eter, X, set equal to zero. 

The APR and Skyrme EOSs are special in that rjpot = 
1 independent of density. In the case of the Skyrme mod- 
els, this results from the structure of the interaction, and 
in the case of APR, r/pot = 1 by construction. 

In Fig. [3 the intrinsic values of 77 (the procedure de- 
scribed in Eqs. I7I9I has not yet been employed) as ob- 
tained from the EOSs described above are plotted, as 
well as that obtained by fitting the results in Ref. 0. 
Note that only the total energy per baryon was reported 
in this latter reference, so there are no results for the top 
panel. The value of T^pot for the typical RMF model is 
quite close to unity, while RMFlo and RMFhi (by con- 
struction) have values significantly different from 1. The 
maximum value of ?7pot for the RMFhi EOS, even though 
it lies outside the graph, is below the limit 7/ < 5 men- 
tioned in the introduction. The MDI EOS has a large 
value of »7potj even at lower densities where most models 
have 77pot ~ 1 

The values of Tytot shown in the bottom panel show a 
little more variation. The value of rytot from the BHF cal- 
culation is quite close to 1, even at lower densities where 
other models (except RMFlo) predict larger values. This 
is typical; in order to get Tytot equal to unity at low den- 
sities, one frequently has 77pot < 1 in order to cancel the 
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effect of Tykin which is most often larger than one. It is 
also quite clear that the effects of the kinetic part of the 
equation of state are not trivial in determining Tytot • The 
ordering of the values of rj for the RMFlo and RMFhi 
EOSs reverses depending on whether or not one is con- 
sidering the total equation of state or only the potential 
energy part. 
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FIG. 2: The density dependence of ?7pot and rjtot determined 
from Eq.|K|for the equations of state considered in this work 
as described in the text. The kink in APR is due to the phase 
transition to the high-density phase. 



V. DIRECT URCA 

When the neutron star temperature is small enough 
that it can be ignored in comparison to the Fermi mo- 
menta of the consituents, the critical density for direct 
Urea is equivalent to the condition that one be able to 
form a triangle with the neutron, proton, and electron 
Fermi momenta. This triangle may be formed if the 
squared area obtained by Heron's formula is positive, i.e. 

s{s - kFn){s - kFp)is ~ kpe) > (11) 

where s is the semiperimeter defined by s = {kpn + 
kpp + kpe)/2. This reduces the familiar condition kpn < 
kpp + kp^ for neutron- rich matter and is more easily gen- 
eralizable to other Urca-like processes. 

When an EOS has a proton fraction which is small, 
true of the APR EOS and some EOSs based on Skyrme 
interactions, the presence of a quartic term can drasti- 
cally affect the critical density for the direct Urea process. 
This is demonstrated in Fig. |31 where the critical density 
is plotted as a function of rjpot , assuming that 77pot is den- 
sity independent. Depending on the relative strength of 
the quartic dependence. Urea either proceeds in all neu- 
tron stars with masses larger than about 1.4 solar masses 



or in no neutron stars of any mass. Although values of 
77pot less than 1/2 are allowed, the corresponding criti- 
cal densities are already above the central density for the 
maximum mass configuration. The maximum mass and 
radius of the maximum mass neutron star as a function 
of tjpot are also plotted in Fig. |21 These quantities are 
essentially unaffected by the modification of ?7pot- 

RMF models typically have a large symmetry energy 
and a large proton fraction, and thus allow the direct 
Urea process at low densities. However, this is not re- 
quired, and Ref. constructed a couple RMF models 
with small neutron star radii (named SRI, SR2, and SR3) 
which do not allow direct Urea at any density lower than 
the central density of the maximum mass configuration. 
These models also typically have a smaller symmetry en- 
ergy at saturation density than the models utilized here. 
The models RMF, RMFlo, and RMFhi, have critical den- 
sities of 0.306, 0.275, and 0.806 fm'^, respectively. Note 
that RMFhi has a rather large critical density for direct 
Urea which could easily be made larger by decreasing the 
value of the symmetry energy at saturation density. 
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FIG. 3; The critical density for the direct Urea process for 
the APR EOS as a function of r^pot- Points with r^pot < 1/2 
were not plotted because the Urea process is not allowed at 
any density for this range of rj. 

Ref. (3J|] compared the radii of 1.4 Mq neutron stars 
with the MDI EOS using different symmetry energies and 
comparing with isospin diffusion data from intermediate- 
energy heavy-ion collisions. This reference found that the 
threshold for the direct URCA process was signficantly 
different for the APR EOS and for the MDI EOS (with 
X = 0) even though they had nearly identical symmetry 
energies. This is, in large part, due to the signficant pres- 
ence of quartic terms in the MDI EOS as demonstrated 
above. This does not contradict the constraint on neu- 
tron star radii from Ref. j34j , since neutron star radii are 
insensitive to quartic terms, as demonstrated in Fig. |51 
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VI. DISCUSSION 

Quartic terms play an important role in determining 
the critical density for the direct Urea process in neutron 
stars. These terms can be easily generated within the 
context of RMF models of high-density nucleonic mat- 
ter. While this work means that it will be more difficult 
to interpret neutron star cooling data without more infor- 
mation on the value of rj at large densities, it also means 
that more neutron star cooling data is essential to under- 
standing the nature of the high-density equation of state. 
Observations of neutron stars masses and radii will have 
difficul ty c onstraining the value of rj. 

Ref. studied in detail the cooling of neutron stars 
constructed with the APR EOS. They found that, be- 
cause of the direct Urea process, stars with masses larger 
than about 1.7 Mq cool sufficiently rapidly as to be cooler 
than nearly all of the observed neutron stars. As these 
authors point out, this is somewhat sensitive to the as- 
sumptions about the pairing interaction, and sufficiently 
cool neutron stars with large mass may be difficult to 
obeserve. Nevertheless, this work offers another possi- 
ble interpretation: the actual value of Typot niay be suffi- 
ciently small at high density to turn off the direct Urea 
process thus making the computed cooling curves match 
the comparatively warm neutron stars which are given in 
the experimental data. 

Proto-neutron star cooling (during the first minute 
when the core is still hot) always contains a significant 
contribution from direct Urea (e.g. see the review in 
Ref. [S^) since the finite temperature allows momentum 
conservation to be easily fulfilled. However, the proton 
fraction of the hot neutron star matter dictates the flavor 
content of the associated neutrino signal and this will be 
modified by the presence of quartic terms. 

Ref. pointed out that there is a correlation between 



the threshold for direct Urea and the neutron skin thick- 
ness of lead. This work clearly indicates that the correla- 
tion is predicated on the implicit assumption that T^pot is 
nearly unity for the relevant densities. This will also, for 
some models, impact the suggested correlation between 
the direct Urea process and neutron star radii 37], since 
quartic terms are important for the former but not for 
the latter. 

It might be interesting to explore the effect of a quartic 
term on the symmetry energy at finite temperature. The 
value of 1] for the kinetic part of the equation of state 
could become closer to unity since finite-temperature ef- 
fects are proportional to T^. Nevertheless, various exper- 
iments are now probing the temperature-dependence of 
the symmetry energy and might provide constraints on a 
quartic behavior [33, • 

The effects of a phase transition to hyperons, Bose- 
condensates (except for that already present in APR), or 
quarks have been ignored. It might be interesting to ex- 
plore the value of r] for equations of state involving phase 
transitions. Although quark matter is often more isospin- 
symmetric, it also frequently has a significant gap, which 
prohibits the direct Urea process for low enough temper- 
atures. 
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